Intermittency induced long-range cross-correlations 
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Cross-correlations are usually considered to emerge through interaction between particles. Here 
we present a mechanism capable to generate power-law cross-correlations between non-interacting 
particles exposed to an external potential. This phenomenon can occur as an ensemble property 
when the external potential induces intermittent dynamics, providing laminar and stochastic phases 
of motion. We have strong indications that the divergence of the mean residence time in the laminar 
phase of the single particle motion - sporadic dynamics - is a sufficient condition for the emergence 
of long-range cross-correlations. We argue that this mechanism may be relevant for the occurrence 
of collective behaviour in critical systems. 
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Complex systems usually consist of several dynamical components interacting in a non-linear fashion. Cross- 
correlations are then used in order to explore the inter-dependence in the time evolution of these components measured 
in terms of specific quantities characterizing each component. In this context, the existence of cross-correlations has 
been demonstrated in a wide class of dynamical systems ranging from nano-devices [I] to atmospheric geophysics 
[2] , seismology [3] , finance |4- 8 , physiology and genomics [8] . Of special interest is the case of long-range (power- 
law) cross-correlations (LRCC) which, being scale free, may be associated with the appearance of characteristics of 
criticality in the dynamics of the considered complex system. Such a behaviour has been observed, among other 
examples, in price fluctuations of the New York Stock Exchange during crisis jS], physiological timeseries of the 
Physiology Sleep Heart Health Study (SHHS) database [8], the spatial sequence describing binding probability of 
DNA-binding proteins to genes at different locations on mouse chromosome 2 [5] and in flocks of birds [5] ■ All these 
findings indicate that the presence of power-law cross-correlations is a quite general property of the dynamics of 
complex systems. Even more, very recently geometry induced power-law cross-correlations have been also observed 
in a coarse-grained description of the dynamics of an ensemble of non-interacting particles propagating in a Lorentz 
channel |10j . This clearly pauses the question of the origin and mechanisms of cross-correlations in particle systems. 

Up to now the theoretical treatment of cross-correlations is based on statistical approaches and their microscopic 
origin is to a large extent unclear. In the present Letter we identify the dynamical mechanisms leading to LRCC 
and show specifically that intermittent dynamics obeyed by each component separately generates LRCC between the 
different components even if they do not interact with each other. In addition we provide strong evidence that a 
sufficient condition for the emergence of such scale-free LRCC is the divergence of the mean length of the laminar 
phase in the intermittent dynamics of each component. 

The prototype model we will use to demonstrate our arguments is a system of N non-interacting particles each with 
an one-dimensional phase space determined by the variable x^ (i — 1, 2, ...,N). We do not further specify iW; in a 
real system it can be for example the position or the momentum of particle i or any other property characterizing its 
state (partially or totally) . We use a version of the well-known Pomeau-Manneville map of the interval [11] , consisting 
a normal form of intermittent dynamics, to describe the time evolution of x^ of each particle independently: 

(«) _/ 4 i) +« w (4 i) )* w ; ^6(0,^*) m 

Xn+1 ~ { r^zW" ; x® € (*(*>.*,!] [i) 

for i = 1,2, N. In Eq.|l]) «W are positive constants, are characteristic exponents fulfilling zW > 1 and r$ are 
random numbers uniformly distributed in (0, 1]. The quantity x^'* represents the upper boarder of the phase space 
region ((0, x^'*)) within which the evolution of the particle dynamics is laminar. Notice that in Eq. M there is no 
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coupling term between phase space variables of different particles since there is no mutual interaction. This simple 
model is quite general capturing all the basic dynamical ingredients necessary for the development of cross-correlations 
as we will see in the following. To avoid unnecessary complexity we further simplify the model assuming: = u 
and = z for all i. Note that the end of the laminar region x* is not strictly defined. One possible choice, which we 
use in the following, is to fix x* as the pre-image of 1, i.e. as the solution of the equation x* + u(x*) z — 1. A second 

choice is to set it equal to x* — y^r) i_1 being the 3>value for which the non-linear term in Eq. (l| becomes equal in 
magnitude with the linear term. These two values (x* and x*) are very close to each other for almost all values of z 
and our results for the cross-correlations, shown below, do not depend on this choice. Using equation ([IJ) we evolve 
the considered particle system in discrete time. Different particles correspond to different trajectories i.e. trajectories 
starting from a different initial condition. Thus we propagate a set of N trajectories. At each point of a trajectory 
we define an observable quantity A„ = A{xn^) (which can be x { n ] itself). Then the cross-correlation function with 
respect to the observable A is defined as: 

CCa ^= N(N-l)^ ^%U) (2) 

where (...) denotes time averaging while and a^jp are the standard deviations of A^ and respectively. 

A typical characteristic of the intermittent dynamics is that for any trajectory the x-values in the laminar region 
are very close to the diagonal x^ +l = Xn since the increase Ax„ = x^\i — Xn of x n ^ there is very slow. It is 
important to notice that the cross-correlations, calculated using Eq. (§, appear in the phase space variable x while 
the increments Aa; of different trajectories are clearly uncorrelated. This behaviour is illustrated in Figs. la,b where 
we plot CC x (m) (Fig. la) and CC^ x {m) (Fig. lb) for an ensemble of trajectories evolving according to the law given 
in Eq. |TJ) with u = 1 and z — 2.5. 

To illuminate further the origin of the cross-correlations shown in Fig. la we explore their dependence on the value 
of the parameter z. For this task we calculate CC x (m) for various values of z using ensembles of 10 4 trajectories with 
length 10 5 for each case. For z > 2 we find an algebraic decay of CC x {m) with increasing m, having an exponent 
which depends on z, while for z < 2 CC x (m) fluctuates around zero indicating the absence of cross-correlations in 
this case. The behaviour of CC x (m) for z < 2 is demonstrated in Fig. 2 where we plot as a typical example the result 
obtained for z = 3/2. It is worth to mention here that a distinction between the properties of intermittent dynamics 
for z < 2 and z > 2 has been already discussed in [12] where the term sporadicity is introduced for the description of 
the z > 2 case. 

In order to explain the different behaviour of the cross-correlation functions for z > 2 and z < 2 we consider 
the distribution of the laminar phase lengths or as it is often also named the distribution of the waiting times in 
the laminar region. It is well known that this distribution obeys asymptotically (I ^> 1) a power-law of the form 
p(£) ~ £~^t [T3J[I3| where £ is the laminar phase length. For z > 2 the mean laminar length (£) diverges while for 
z < 2 it is finite [T7] . Thus the divergence of (£) should be related with the emergence of cross-correlations between the 
particles. This in fact can be explained as follows: if {£) is infinite then the conditional probability that the particle 
j at an instance n + m is in the laminar region provided that the particle i was in the laminar region at instance n is 
finite and decays algebraically with increasing m. 

To understand better this behaviour it is useful to develop a symbolic code for the intermittent dynamics in 
Eq. ([!]). Such a symbolic representation of the dynamics of the Pomeau-Manneville intermittent map capturing 
several details like the re-injection rate in the laminar region (and therefore the invariant density in the immediate 
neighbourhood of the marginally unstable fixed point) is proposed in [T^]. Here we are interested mainly to isolate the 
dynamical properties leading to the emergence of cross-correlations avoiding the influence of other detailed aspects of 
the intermittent dynamics. Therefore we will use a much simpler code, mapping x in the laminar region (x € [0,a;*]) 
to and x out of the laminar region (x € (a;*,l]) to 1. Such a code is used in 14 to calculate power-spectra 
of intermittent systems. In practice we use the full dynamics of Eq. to generate the ensemble of intermittent 
trajectories and then we replace the x-values in each time-series by or 1 according to the previously described rule. 
Subsequently we calculate the cross-correlation function CC s (m) for different z-values using the binary sequences 
generated by the symbolic dynamics from the trajectories of the map in Eq. ([T]). 

Furthermore, in order to get rid of any other dynamical effects which may still be present in the symbolic dynamics, 
we introduce a stochastic model containing only the information of the laminar length distribution to simulate the 
emergence of cross-correlations. We assume a process consisting of two phases defined as follows: (i) a stochastic 
variable £ takes the value 1 in the irregular phase and the value in the laminar phase and (ii) the length of the 
irregular phase is always 1 while the laminar length probability density is a power-law with exponent —zj(z — 1) (z 
being the exponent in the intermittent map of Eq. 0). Then we generate an ensemble of realizations of this process 
and calculate the cross-correlation function CC' r (m) for this ensemble. Despite the simple form of the map providing 
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FIG. 1: (color online), (a) The cross-correlation function CC x (m) as denned by Eq. ^ using 10 4 trajectories each of length 
10 5 and, (b) CC A x{m) for the same ensemble of trajectories. The parameters of the intermittent dynamics are z = 2.5 and 
u = 1. 
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FIG. 2: (color online). The cross-correlation function CC x {m) for the intermittent map of Eq. |Tj with parameters z = 1.5 
and u = 1. For the numerical simulations we used an ensemble 10 4 trajectories each of length 10 5 . 
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the intermittent dynamics, large scale computational efforts (10 5 trajectories have been propagated for 10 6 iterations) 
are needed to achieve convergence of the long-time behaviour of the cross-correlation function. In Fig. 3 we show the 
results obtained for CCa(w) with A = s,r for z = 2.5, 3, 4, 5. The coloured triangles correspond to A = s while the 
red lines to A = r. We observe a very good agreement between the two results for each z value. This is a strong 
indication that the quantity determining the properties of the cross-correlation function is indeed the laminar length 
distribution. 
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FIG. 3: (color online). The cross-correlation function CCA(rn) for the symbolic dynamics generated from the intermittent 
map of Eq. [A = s, triangles) and for the stochastic process defined in the text (^4 = r, red lines) for four different values 
of z: z — 2.5 (black), z = 3 (blue), z = 4 (green), z — 5 (dark yellow). For the numerical simulations we used an ensemble 10 5 
trajectories each of length 10 6 . 

Based on the stochastic model introduced previously one can develop a method to find an analytical estimation 
CC a (m) of the cross-correlation function CC r {m). Such an attempt will enlighten further the origin of the found cross- 
correlations. Let us first consider two binary sequences {x^} — {xi,x^ \ ■■} and {x^} = {x^jX^ , -,Xn , ■•} 
generated by the stochastic model. The function C'C a (rn) should be proportional to the joint probability Pijfau — 

1; x^)_ m = 1) that the random variable x^ has the value 1 at time step k and the random variable x^ has the value 
1 at time step k + m, averaged over the time: 

1 N—m 

CCa{m) = JT^n £ = = !) ® 

fe=i 

Obviously it holds: 

PiA*f = lisjg™ - 1) - = 1) ' P(4irn = 1) (4) 

since a;W and x^ are statistically independent. To calculate P(x^ = 1) one can use the method introduced in [T4"] 
writing: 

p(4 ] = i) = J2 p ( x k-n = i) • ahh41 = i) (5) 

71=1 

where Pi\i(n\x^_ n = 1) is the conditional probability to have a laminar phase of length n directly after the instant 
fc — n if has the value 1 at the time instant fc — n. The appearance of a laminar phase with duration n is independent 
of the value of at the instant k — n. Thus we find: 

Pi\M4-n = l ) = P( n ) ; P (n)~n-^ ; n»l (6) 
where p(n) is the laminar length distribution. Inserting Eq. ([6| into Eq. ^ we obtain the equation 

fe-i 

P(xt ) = l) = ^2P(4l n = l)-p(n) (7) 

71=1 
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which can be solved recursively using as initial condition P(x{ = 1) = po with po € (0, 1). A similar equation is 
obtained also for P(%^+ m = 1) replacing simply k with k + m. Having solved Eq. |?| one can then calculate the sum in 
Eq. Q to obtain an analytical expression for CC Q (m) containing three sums. The validity of the introduced analytical 
scheme is tested in Fig. 4 where we show the symbolic dynamics result CC s (m) together with the analytical form 
CC a (m) for z = 3. We observe a very good agreement between the analytical result and the numerical simulations 
for CC s (m). Notice that CC r (m) is not shown in this plot. However, as illustrated in Fig. 3 the results for CC r (m) 
and CC s (m) are very close to each other for any considered z and therefore CC a (m) can be used as an analytical 
estimation of the cross-correlation CC s (m) too. The analytical treatment leads us to the conclusion that it is the 

long-range character of the correlation between P{x$ = 1) and P(xff = 1) existing for any pair of intermittent 
trajectories which generates the observed cross-correlations. Note that this property has been discussed in [TS] in a 
different context. 




FIG. 4: (color online). The cross-correlation functions CC s {m) (symbolic dynamics, intermittent map, black circles) and the 
analytical estimation CC a (m) (red line) for z = 3. 

With our analysis we have demonstrated a mechanism to establish power-law cross-correlations between particles 
which do not interact with each other. This phenomenon is induced by the strong intermittent dynamics performed 
by each of the particles independently. Strong intermittency (sporadicity) discussed in this Letter is a result of 
the interaction of a particle with a suitable external potential (field) [TS]. The appearance of long-ranged cross- 
correlations deems sporadic dynamics a plausible mechanism for the collective behaviour emerging in a A-particle 
system. Furthermore, since such a collective behaviour is accompanied by scale free inter-particle correlations, it 
could be related to the emergence of critical behaviour in the considered system. In fact a connection of intermittent 
dynamics with criticality has already been established in [T^] using the example of the 3 — 1? Ising model. There it has 
been shown that the order parameter fluctuations at the critical point can be efficiently described by an intermittent 
map of Pomeau-Manneville type (similar to that of Eq. Q) with additive noise. The exponent z in this intermittent 
map is related to the isothermal critical exponent 8 associated with the second order transition. This property sets a 
bound z > 2 necessary for the occurrence of critical behaviour. It is remarkable that this bound coincides with the 
bound obtained by our present analysis in order to have a divergent mean laminar length. An astonishing feature 
of our results is that the power-law cross-correlations emerge even without interaction among the particles. In the 
context of critical phenomena such a property is welcome since it could explain universality aspects. Indeed the 
microscopic interactions between the elementary degrees of freedom of a critical system do not play any role for the 
determination of the critical exponents and the associated scaling laws describing the phenomenology of an extended 
system at the critical point. 

In the framework of our approach the obtained correlations are determined by the time evolution of the trajectories 
of two different particles. To enable a closer relation to equilibrium critical phenomena one should extend these ideas 
also to the case of a field depending both on time and on space. Such an extension requires the use of matrix equations 
for the field evolution replacing the variable Xn by a scalar field (f)(i,n) where i is a spatial variable while n is the 
time variable. At a first glance one could argue that for the calculation of the spatial cross-correlations one might 
exchange the role of spatial and temporal variables in the dynamics, use Eq. ([I]) to describe changes of the field <f> in 
space and average over the time variable. This would lead to power-law cross-correlations between the field values at 
different locations which is typical for a critical system. However a consistent treatment of this case requires more 
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elaborate and extensive studies left for future investigations. 
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